Variation of Specific Heat with x and y in Na x Co02 ?/H20/D20 
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We report specific heat measurements down to 0.4 K on the layered oxide Na^CoC^ • 2/H2O/D2O 
with < x < 0.74 and y — and 1.4. For the nonhydrated system (y — 0), the electronic 
specific heat coefficient 7jv and the Debye temperature 6d vary nonmonotonically with x, both 
displaying minima when x is close to 0.5. This indicates a systematic change of the electronic and 
vibrational structures with Na content. For both hydrated and deuterated systems (x = 0.35 and 
y = 1.4), the specific heat reveals a sharp peak with AC P ~ 45.5 mJ/mol-K at T™ ~ 4.7 K and 
an anomaly at T x ~ 0.8 K. While the origin of the later is unknown, the former corresponds to the 
superconducting transition. Interestingly, the electronic specific heat, after subtracting lattice and 
Schottky contributions, exhibits roughly T 2 behavior between 0.2T C and T c . This can be explained 
by assuming an unconventional superconducting symmetry with line nodes. The results obtained 
under applied magnetic field further support this scenario. 



PACS numbers: 74.20.Rp, 74.25. Bt, 74.25.Jb, 74.90. +n 



There is growing evidence that the strong electron- 
electron correlation in layered Na^CoC^ is responsible 
for some of its anomalous physical properties such as its 
"colossal" thermopower 0,0 and possibly superconduc- 
tivity in its hydrated form.Q Although it is a good metal 
with high electrical conductivity for a wide Na doping 
range except for x — 0.5, both local density approxima- 
tion (LDA) calculations and experimental work [j| 
indicate that the itinerant bands of Na^OoC^ are very 
narrow with W <C U, where W is the band width and 
U is the effective on-site Coulomb interaction. This im- 
plies a high value of the density of states (DOS) at the 
Fermi level. However, all specific heat data reported so 
far reveal a weak or moderate enhancement of electronic 
specific heat coefficient 77V for both hydrated and non- 
hydrated cases compared to the value from LDA band 
structure. 0, El IS 01 In these reports, the 7 at value was 
obtained by analyzing specific heat data above 2 K for 
Na content x in the range of 0.3 - 0.8. It is possible that 
the extracted value does not represent that for T w 
K. On the other hand, recent calculations, using the 
LDA+U method, suggest that the strength and effect of 
Coulomb interactions are reduced with decreasing x.[Toj 
One would thus expect the variation of 7 at with x. 

On the experimental side, the electronic properties of 
hydrated and nonhydrated Naa;Co02 have not been sys- 
tematically studied as a function of composition. While 
the phase diagram shown in Ref. is constructed from 
electrical transport and magnetic measurements for 1/3 
< x < 3/4 and y = 1.4, little is known about how the 
thermodynamic properties vary with both x and y. Of 
particular importance is the specific heat behavior in the 
superconducting state of the system, which can provide 
key information about the superconducting pairing sym- 
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metry. In this Letter, we report the low-temperature 
specific heat of Na K Co0 2 • j/H 2 0/D 2 with < x < 0.74 
and y = and 1.4. 

Single crystals of Na x Co02 were used for specific heat 
measurements. Starting with Nao.74Co02 single crystals 
grown using a flux method. |l2j crystals with lower Na 
content were obtained by chemical dcintcrcalation as de- 
scribed in Ref. [fj. By controlling the deintercalation 
time, we obtain single crystals with x ss 0.72 and 0.48, 
as determined from measurements of the c-axis lattice 
parameter using the calibration curve in Ref. [llj . How- 
ever, it is known that the above technique cannot ex- 
tract all Na from Naa;Co02. In order to obtain C0O2 
(.x = 0) with the hexagonal structure, we extract all 
Li from LiCo02 powders using N0 2 BF4 as described in 
Ref. [LJ- Polycrystalline samples of Nao. 35 C0O2 and su- 
perconducting Nao.35CoO2-l.4H2O/D2O were prepared 
following a procedure similar to that described in Ref. 
|14|. Specific heat measurements were carried out using 
a PPMS (Physical Property Measurement System) from 
Quantum Design. 

In Fig. la, we present the temperature dependence 
of the specific heat C p for single crystalline Nao.74Co02 
(filled circles), Nao.72Co02 (unfilled circles), Nao.4sCo02 
(crosses), polycrystalline Nao.3sCo02 (filled diamonds), 
and C0O2 (unfilled diamonds). Note that C p varies 
monotonically with T between 0.4 and 10 K and the curve 
tends to move upward with increasing x. For easy analy- 
sis, we replot the data as C p /T versus T 2 as shown in Fig. 
lb. If the specific heat is due to electrons and phonons 
only, it is expected that C p /T will be proportional to 
T 2 at low temperatures. While this is true for C0O2 (x 
~ 0) between 0.4 and 10 K, C p /T clearly deviates from 
linearity below ~ 2 K for samples with x ^ 0. The low- 
temperature upturn indicates that an additional contri- 
bution sets in which increases with decreasing T. Similar 
behavior was reported previously and was attributed to 
a Schottky effect. Thus, for T < Q D (Debye temper- 



2 




T(K) 7 s (K 2 ) 

FIG. 1: (a) Temperature dependence of the specific 
heat C p for single crystal Nao.74Co02 (filled circles), 
Nao.72Co02 (unfilled circles) , Nao.4sCo02 (crosses), polycrys- 
talline Nao.3sCo02 (filled diamonds), and C0O2 (unfilled di- 
amonds) between 0.4 and 10 K. Shown in (b) is the replot of 
the data as C p /T versus T 2 , and solid curves are the fits of 
experimental data to Eq. 1. 



ature), we may describe the specific heat using 

C p /T = lN + (p 3 T 2 + f3 5 T 4 )+a £T (1) 

1 J (l + eT y 

Note these terms are corresponding to contributions 
from electrons, phonons (terms inside the bracket), 
and Schottky anomalies, [15| respectively. Here, (3-$ = 
N(12/5)7r 4 Re^ 3 with R = 8.314 J/mol-K and N = 
atomic number per formula unit, f3§ is a constant de- 
scribing the anharmonic coupling, a is a constant pro- 
portional to number of two-level systems, and A is the 
energy separation between the two levels. 

The solid curves shown in Fig. lb are fits of experi- 
mental data between 0.4 and 10 K to Eq. 1. Note that 
Eq. 1 describes our experimental data very well in the 
selected fitting range for all studied compounds. Param- 
eters obtained from the above fitting procedure are listed 
in Table I. Using the relationship described above, the 
Debye temperatures may be calculated from the ^3 val- 
ues, and are given in Table I. In view of Table I, there 
are a couple of remarkable features: (1) both and 
03 (6d) vary nonmonotonically with x, displaying ex- 
trema for x ~ 0.48; (2) a decreases with decreasing x; 
(3) A does not vary much with x. Compared to the cal- 
culated band value of 7^ = 10 mJ/mol-K 2 for x > 0.5 
HE3 and 13 mJ/mol-K 2 for x < 0.50, our 7jv val- 
ues have a factor of 3 enhancement for x > 0.48 but are 
very close for x < 0.48. This implies that the correlation 
effects are moderate for x > 0.48 and become negligi- 
ble for x < 0.48, in very good agreement with LDA+U 
calculations. 0| The small 7 at value for x ~ 0.48 is likely 
related to the charge ordering reported in Nao.sCoC^ at 
temperatures below 52 K. 11 For Nao.4sCo02, there may 
exist partial charge ordering, leading to a small 77V value. 
It should also be mentioned that the value for the end 
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TABLE I: Parameters obtained from the fit of experimental 
data between 0.4 and 10 K to Eq. 1 and the inferred Debye 
temperatures Qd (see text) for all x values. 



compound C0O2 is nonzero, consistent with the metallic 
behavior reported previously. Interestingly, the vari- 
ation of Na concentration in Na^OoC^ affects not only 
the DOS (reflected by 7^) but also @£> with a minimum 
occurring at the same composition (x ~ 0.48) as that 
for 7tv. This suggests strong coupling between electrons 
and the lattice, and the lattice tends to soften near x = 
0.5. The monotonic decrease of a with x reflects that 
the Schottky contribution to the specific heat decreases 
with Na content, and is almost negligible in the end com- 
pound C0O2. This is consistent with the trend that the 
low-temperature magnetic susceptibility decreases with 
decreasing The almost constant A value for dif- 

ferent x strongly suggests that Schottky effect is mainly 
controlled by nuclear moments. 

We now focus on the specific heat of hy- 
drated/deuturated Na .35CoO2. Shown in Fig. 
2a is the temperature dependence of specific 
heat for Nao.35CoO2-l.4H2O (filled circles) and 
Nao.35CoO2T.4D2O (unfilled circles) between 0.4 and 
10 K. Compared to that of the parent compound (filled 
diamonds), the specific heat of hydrated/deturated 
Nao.35Co02 is larger at high temperatures but smaller 
below ~ 2 K. Most prominent is the specific heat 
anomaly below 5 K in both hydrated and deuturated 
cases, reflecting the superconducting phase transition 
at T° nset ~ 5 K. To obtain superconducting-state prop- 
erties, all non-electronic contributions, which are not 
affected by the superconducting phase transition, should 
be subtracted from the total specific heat. As shown in 
Fig. 2b, Cp/T reveals more or less linear dependence 
with T 2 above T c for both hydrated and deuturated 
compounds. This suggests that the first two terms of Eq. 
1 are dominating. While the low-temperature upturn is 
absent, the contribution from the Schottky form may be 
hidden due to sharp decrease of electronic specific heat 
in superconducting state. For comparison, we extract 
the electronic specific heat with/ without considering 
the Schottky effect. Shown in Figs. 2c-2d are the 
temperature dependence of the electronic specific heat 
assuming Cf = Cp-Cf (Fig. 2c) or C e p l = Cp-Cf-G^ 11 
(Fig. 2d). Here, parameters for are assumed to 

be the same as that for Nao.35Co02 (see Table 1) since 
they cannot be determined by fitting data above T c , 
and C|J (T) are obtained by fitting experimental data 



3 



0.4 



o 

E0.2 



0.0 



• Na .35CoO 2 -1 .4H 2 a° 
o Nao.35CoO2l.4D2O . 

♦ Nao. 35 Co02 



(a 




4 6 8 10 20 40 60 80 100 
r(K) T 2 (K 2 ) 





25 
20 
15 

H 

c\j 
u 

a: 10 

5 




• Nao.35CoO2l.4H2O 
o Nao.35CoO2l.4D2O 
■■ WHH 



• Nao.35CoO2l.4H2O 
o Nao.35CoO2l.4D2O 

o 



1 



2 3 
T(K) 



4 5 




1 2 3 
H°-5 (T° ^ 



FIG. 2: (a) Temperature dependence of the specific heat 
C p for polycrystalline Nao.35CoO2T.4H2O (filled circles), 
Nao.35CoO2T.4D2O (unfilled circles), and Na .35CoO 2 (filled 
diamonds) between 0.4 and 10 K; (b) replot of the data in 
(a) as Cp/T versus T 2 ; solid curves are fits of experimental 
data to Eq. 1 (see text); (c) and (d) present the temperature 
dependence of (C p -Cf)/T (c) and (C p -C p p h -C^ ch )/T (d) of 
Nao.35CoO2T.4H2O (filled circles) and Nao.35CoO2T.4D2O 
(unfilled circles). In each case, the solid line is the linear fit 
of experimental data between 0.8 and 4.1 K and the broken 
line is the extrapolation of experimental data below 0.7 K to 
the origin. 



FIG. 3: (a) Temperature dependence of the specific heat 
C p for polycrystalline Nao.35CoO2T.4H2O under indicated 
magnetic fields; (b) temperature dependence of the elec- 
tronic specific heat with Cp /T=(Cp-C^-Cp C,l )/T; (c) tem- 
perature dependence of the upper critical field H C 2 for 
Nao.35CoO2T.4H2O (filled circles) and Nao.35CoO2T.4D2O 
(unfilled circles). The dashed line is the fit of the data to 
the WHH formula (see text); (d) magnetic field dependence 
of the residual electronic specific heat coefficient 7 re s plotted 
as j res vs. H ' 5 . The solid lines are the linear fit of the data 
(see text). 



between 6 and 10 K without (Fig. 2c) or with (Fig. 2d) 
the consideration of the Schottky effect. It turns out that 
the results are almost the same at high temperatures 
with only a slight difference at low temperatures. In 
both cases, note that the substitution of D for H has 
little effect on the specific heat as we concluded earlier 
in Ref.Jfj. As marked by the dashed line in Figs. 2c-2d, 
the normal-state electronic specific heat coefficient 
77V=13.7 mJ/mol-K 2 for Nao.35CoO2T.4H2O and 13.3 
mJ/mol-K 2 for Nao.35CoO2T.4D2O, is very close to that 
for the parent compound Nao.35Co02 (see Table 1). This 
is in very good agreement with a theoretical prediction 
that the intercalation of water into Nao.35Co02 results 
in little change in the electronic structure, [lq 

Although the specific heat starts to depart from the 
high-temperature behavior at 5 K and peaks at ~ 
4.1 K, the thermodynamic transition temperature of 
Nao.35CoO2T.4H2O/D2O is at T™ d - 4.7 K where spe- 
cific heat jump reaches half of the maximum. By extrap- 



olating the low-temperature specific heat to T™ lrf (see 
Figs. 2c-2d), we obtain the specific heat jump AC P ~ 45.5 
mJ/mol-K. It follows that AC p /j N T]! lid ~ 0.71 mJ/mol- 
K 2 , roughly 50% of the BCS value. According to previous 
work, this result is compatible with (a) 50% supercon- 
ducting volume fraction [H [3 or (b) 100% supercon- 
ductivity if there exists line nodes in the superconducting 
pairing symmetry. [TtI IT9I ] As the specific heat measure- 
ments are always carried out in high vacuum, it is possi- 
ble that some H2O/D2O is pumped out of the material 
thus reducing the superconducting volume fraction. For 
this reason, samples were cooled rapidly to 20 K in He- 
lium gas before pumping. Thermogravimetric analysis 
indicates that the H2O/D2O content y remains 1.4 after 
the specific heat measurements. This and the fact that 
all othe r sp ecific heat measurements on either powder 
[H El Ea or single crystals @] result in AC p /7 W T™ d 
up to 50% of the BCS value indicate that the small spe- 
cific heat jump is unlikely due to the sample quality. In 
the second scenario, one would expect a T 2 dependence of 
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the electronic s peci fic heat in the superconducting state, 
i.e. C;'/T oc T. 20] As may be seen in Figs. 2c-2d, C e p l /T 
varies more or less linearly between 0.8 and 4.1 K. The 
solid lines are the linear fit to the experimental data in 
this region. A somewhat steeper slope is obtained when 
the Schottky effect is considered (see Fig. 2d). Thus, 
it is conceivable that there exist line nodes in supercon- 
ducting pairing symmetry of Nao.35CoO2T.4H2O/D2O. 
However, /T shows an apparent deviation from linear- 
ity with sharper decrease below T x ~ 0.8 K. This implies 
that not all of the conducting electrons are involved in 
the superconducting phase transition at T c , leading to 
a small specific heat jump and a non-zero intercept j res 
when the solid lines in Fig. 2c-2d are extrapolated to T 
= K. Residual unpaired electrons should be responsible 
for the anomaly at T x , below which C p /T seems to ap- 
proach linearly to the origin as guided by the broken lines. 
It should be mentioned that the entropy is well balanced 
when the Schottky contribution is taken into account (see 
Fig. 2d). This strongly suggests that the sharp decrease 
of electronic specific heat below T x is unlikely due to ex- 
perimental error. I2H Although the specific heat reported 
by Oeschler et al [ 231 was measured above 0.8 K, it was 
argued in Ref. :22] that the temperature dependence of 
specific heat of Nao.35CoO2T.4H2O should behave sim- 
ilarly to that of MgB2 with two-gap superconductivity, 
where a sharp decrease of specific heat well below T c is 
observed due to a small-gap contribution. 

To gain insight into the transition at T c and 
anomaly at T x , we have measured the specific heat of 
Nao.35CoO2T.4H2O and Nao.35CoO2T.4D2O under the 
application of magnetic field. Shown in Fig. 3a is spe- 
cific heat data on Nao.35CoO2T.4H2O plotted as C p /T 
vs. T with H = (red), 1 T (pink), 2 T (orange), 4 
T (forest green), 7 T (green), 10 T (blue), and 14 T 
(black). Interestingly, the anomaly at T^ seems to be 
insensitive to H. This suggests that the anomaly at T x is 
unlikely due to the small-gap contribution in the two-gap 
superconductivity scenario. |22j As demonstrated in Fig. 
3b, the anomaly remains after subtracting the phonon 
and Schottky contributions. To determine whether it 
is related to two-gap superconductivity or some sort of 
magnetic transition such as a spin density wave, further 
investigation is necessary. 

Nevertheless, the specific heat jump due to sper- 
conducting phase transition weakens and moves to- 
ward lower temperatures with increasing H, indicat- 
ing the suppression of the superconducting transi- 
tion temperature and superconducting volume frac- 
tion. In Fig. 3c, we plot the upper critical field H C 2 



versus T for Nao.35CoO2T.4H2O (filled circles) and 
Nao.35CoO2T.4D2O (unfilled circles). Compared to that 
obtained from the electrical resistivity^ and magnetic 
susceptibility measurements . |23j H C 2(T) shown in Fig. 
3c varies moderately with initial slope dH C 2/dT|r c ~ - 
5.3 T/K. Using this value and the Werthamer-Helfand- 
Hohenberg (WHH) formula for the dirty limit, [24} we 
may estimate the temperature dependence of H C 2 be- 
tween and T c . Remarkably, our experimental H C 2(T) 
follows very well the WHH expression represented by the 
dashed line in Fig. 3c. At T = K, the WHH formula 
gives H C 2(0) = 17.1 T. This allows us to estimate the co- 
herence length £(0) ~ 44 A using H C 2 = 4>o/2-k£ 2 (6o is 
the flux quantum). As emphasized previously. |3l l23| the 
coherence length is short for a superconductor with such 
a low T c . On the other hand, the reduction of the super- 
conducting volume fraction due to the application of H 
is reflected by the increase of 7 res vaule. From linear ex- 
trapolation of Cp /T(T) in the linear regime of supercon- 
ducting state, we obtain 7 res for each H and plot them 
as 7 res versus H 05 in Fig. 3d for Nao.35CoO2T.4H2O 
(filled circles) and Nao.35CoO2T.4D2O (unfilled circles). 
As guided by the solid lines, 7 res varies more or less lin- 
early with H 0,5 for both systems, consistent with uncon- 
ventional superconducting pairing symmetry with line 
nodes. 

Based on the present specific heat studies, we conclude 
that the electronic and vibrational properties are sensi- 
tive to x but not y in Naa,Co02 ■ 2/H2O/D2O. For the 
nonhydrated system with y = 0, both the electronic spe- 
cific heat coefficient 7jv and the Debye temperature Qd 
change systematically with x, showing minima at 1 ~ 
0.48. For the hydrated/deturated system with x = 0.35, 
the variation of y between and 1.4 results in little effect 
on 7jv value. There is evidence for "double transitions" 
with one at T c ~ 4.7 K and another at T x ~ 0.8 K. While 
the origin of the later "transition" is unknown, our re- 
sults demonstrate that the superconducting transition at 
T c is likely unconventional with line nodes in its pairing 
symmetry. 
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